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Abstract 

The nuclear symmetry energy coefficient (including the coefficient OsyL of D term) of finite nuclei 
is extracted by using the differences of available experimental binding energies of isobaric nuclei. 
It is found that the extracted symmetry energy coefficient aly^{A,I) decreases with increasing of 
isospin asymmetry I, which is mainly caused by Wigner correction, since e*yj^ is the summation 
of the traditional symmetry energy Csym and the Wigner energy ew- We obtain the optimal values 
J = 30.25 ±0.10 MeV, Uss = 56.18 ± 1.25 MeV, aiyL = 8.33 ± 1.21 MeV and the Wigner parameter 
X = 2.38 ± 0.12 through the polynomial fit to 2240 measured binding energies for nuclei with 
20 < ^ < 261 with an rms deviation of 23.42 keV. We also find that the volume symmetry 
coefficient J ~ 30 MeV is insensitive to the value x, whereas the surface symmetry coefficient Oss 
and the coefficient aiym are very sensitive to the value of x in the range 1 < x < 4. The contribution 
of aiyln term increases rapidly with increasing of isospin asymmetry I. For very neutron-rich nuclei, 
the contribution of aiyL term will play an important role. 
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I. INTRODUCTION 


It is evident that the symmetry energy coefficient plays an extremely important role, not 
only in nuclear physics, such as the dynamics of heavy-ion collisions induced by radioactive 
beams, the proper description of the nuclear binding energies along the periodic table, and 
the structure of exotic nuclei near the nuclear drip lines [1-8], but also a number of important 
issues in astrophysics, such as the dynamical evolution of the core collapse of a massive star 
and the associated explosive nucleosynthesis [9-15]. In the global htting of the nuclear 
masses in the framework of the liquid-drop mass formula, the symmetry energy per particle 
is usually written as Csym = CLsymP, in which the symmetry energy coefficient Osym enters 
as a mass-dependent phenomenological parameter [16-20]. In fact, Ogym could also be a 
function of the isospin asymmetry I = [N — Z)/A. The isospin dependence of the symmetry 
coefficient Osym is usually written as asym(^,-f) = J — + oiyL/^ by neglecting the 

higher order term, the same as in Ref. [21, 22]. But how to change it depends on the isospin 
asymmetry I for given mass number A, decreases or increases? It is mainly determined by 
the high-order term coefficient asyL of the symmetry energy. However the coefficient OsyL 
is difficult to be determined. It is necessary to investigate the symmetry energy coefficient 
of hnite nuclei. 

In Ref. [23], Min Liu et ah obtained the mass dependence of asym(A) through performing 
a two-parameter parabola htting to the energy per particles after removing the Coulomb 
energy e„(H,/) = e(A,I) — ec(A,I) for a series of nuclei with the same mass number A. 
The extracted Usym is only dependent on mass number A. In this work, with the similar 
approach in Ref. [23], we consider the mass and isospin dependence of Osym, and at the 
same time the higher-order (/^) term of the symmetry energy is included. It is found that 
the Wigner energy Eiy should be considered in the extraction of nuclear symmetry energy 
coefficient. The Wigner energy can be extracted from the difference of e„(H, J) of isobaric 
nuclei. However the Wigner energy is not included in extracting the symmetry energy 
coefficient in our previous paper [24]. The nature of the symmetry and Wigner energy are 
intertwined in the nuclear mass formula and that one term cannot be reliably determined 
without knowledge of the other [25]. This leads to considerable uncertainty in the value for 
the symmetry energy, especially the coefficient aiyL of the term in the symmetry energy 
coefficient expression. 
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The paper is organized as follows. In Sec. II, the symmetry energy and Wigner energy 
are described and the summation of both are extracted by using the experimental binding 
energies differences between isobaric nuclei. In Sec. Ill, The method of extracting the 
symmetry energy coefficient is described and the corresponding coefficients are obtained 
through the polynomial htting. The effect of the Wigner energy term on the symmetry 
energy coefficient is also studied, in Sec. IV. Finally a summary is given in Sec. V. 


II. SYMMETRY ENERGY AND WIGNER ENERGY 

It is well known that nuclear mass is one of the most precise measured quantity in nuclear 
physics. It can provide information of the symmetry energy coefficient through the liquid- 
drop mass systematics. In semi-empirical Bethe-Weizsacker mass formula [26, 27], the energy 
per particle e(A, I) of a nucleus can be expressed as a function of mass number A and isospin 
asymmetry I = {N — Z)/A, 

1) = O'v ctsA -\r ec(A, /) -|- asyml'^ + S-, (1) 

with 

6 = ±apA~^^‘^ or 0, (2) 

where the “-I-” is for even-even nuclides, the is for odd-odd nuclides, and for odd-A 
nuclides (i.e. even-odd and odd-even) <5 = 0. The a^, a^, agym and a^ are the volume, surface, 
symmetry and pairing energy coefficients, respectively. The Coulomb energy per particle is 
ec{A, I) = Ec/A, where the Coulomb energy of a nucleus = 0.71-^^(1 — 0.76Z^/^) and 
Z = ^{1 — I) are usually used [28, 29]. 

Let us assume the binding energy per particle e{A,I) = en{A,I) -f- ec{A,I), en{A,I) 
and ec(A, I) denote the nuclear energy part and the Coulomb energy part per particle, 
respectively. Subtracting the Coulomb energy term from the binding energy, one obtains 
the nuclear energy part per particle, 

en{A,I) = e(A,J)-e,(A,J) 

= eo(A)-7 esym(A,/) 

= ^o{A) + asyra{A, I)P, (3) 
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where Cq^A) = + asA~^/^ + 5 including the volume, surface and pairing energy terms, is 

only dependent on nuclear mass number A. esyra{A, I) is the symmetry energy per particle of 
a nucleus. If we take the difference in the nuclear energy part per particle en{A, I) between 
two isobaric nuclei with same odd-even parity, eo(A) term is canceled and the difference of 
the symmetry energy per particle can be written as 

Ae^ym = en{A, I) - e„(A, h) = a^ymiA, J)/^ - a^ym{A, (4) 

Here en{A, Ji) is the nuclear energy part per particle of a reference nucleus (H, Ji), and the 
symmetric nuclei (Ji = 0) is selected as the reference nucleus if its experimental binding 
energy is exist for even-even nuclei. For any other case, the nuclei with the minimum value 
of Ji = Imin > 0 is selected as the reference nucleus among each series isobaric nuclei. 
en{A, I) is the any other value of isobaric nuclei for given mass number A. 

If the experimental binding energy of a symmetric nucleus (/i=0) is known, we obtain 

esym(-4, 1) = e„(H, I) - en(H, 0) = a^ymiA, J)/^ (5) 

or 

asym(4l, I) = - -p - = --, (6) 

where only even-even nuclei are taken into account in our calculations to consider the paring 
effects for the even mass number nuclei. 

On the other hand, according to the liquid drop model, the symmetry energy coefficient 
of a hnite nucleus is usually written as 

asym(^, I) = a^yli + 

= ~ J - assA-^^^ + (7) 

by using the Leptodermous expansion in terms of powers of J 28 —34 MeV denotes 

the symmetry energy of nuclear matter at normal density, a^s is the coefficient of the surface 
symmetry term, aiyh is the coefficient of the term in the expression of symmetry energy. 

Figure 1 shows the extracted experimental symmetry energy coefficients as a function 
of isospin asymmetry I extracted from Eq. (6) for all even-even nuclei with mass number 
H=80 (solid squares), where en(A,I) = e(A,I) — ec(A,I), the experimental binding energy 
per particle e(A, I) is taken from the mass table AME2012 [30], and ec(A, I) = 0.71^p^(l — 
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FIG. 1: (Color online)Experimental symmetry energy coefficients as a function of I extracted 
from Eq. (6) for all even-even nuclei with mass number ^=80 (solid squares). The dotted line 
(asym=0), the solid line (aiym=50 MeV) and the dashed line (aiym=-50 MeV) are the results using 
the expression of symmetry energy coefficient of Eq. (7). 

0.76Z^/^). The dotted line (aiyln=0), the solid line (aiyln=50 MeV) and the dashed line 
(osym = —50 MeV) are the results using the expression of symmetry energy coefficient Eq. 

/o\ 

(7) with asym=23 MeV. From hgure 1, one can see that only using the expression of symmetry 
energy coefficient of Eq. (7), the extracted experimental symmetry-energy coefficient can 
not be reproduced whatever it is positive, zero or negative value for aiyL- 

The effect of the Wigner energy is responsible for the decrease of esym{A,I)/P with 
isospin asymmetry / at a given mass number A. To reproduce the experimental data better, 
one should include the Wigner energy term in Eq. (5). Let us rewrite the expression of 
Eq. (5) as = en{A,I) — en{A,0), where e*y^{A,I) is dehned as the summation 

of the traditional symmetry energy esym(^,/) and the Wigner energy ew{A,I). However 
the different Wigner energy expression and parameters will directly affect the extraction of 
symmetry energy coefficients. Figure 2 (a) presents two forms for Wigner energy which is a 
function of isospin asymmetry I and applied to all even-even nuclei with mass number H=80 
in mass table AME2012. One is ew = 29.156/^[(2 — \I\)/{2 -|- |/|A)] (solid triangles), which 
is proposed in Ref. [28], the other is ew = —10exp(—4.2|J|)/H [32](sohd circles), which is 
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usually used in the literature. For convenience we denote the former by “form (1)” and the 
latter by “form (2)”, respectively. From Fig. 2 (a) one can see that the value of ew is positive 
for form (1) and negative for form (2). While the value is the summation of the 

traditional symmetry energy and the Wigner energy, negative Wigner energy of form (2) will 
lead to a larger traditional symmetry energy and thus larger symmetry energy coefficient 
than that with form (1). Fig. 2 (b) presents the extracted symmetry-energy coefficients 
Osym by using two Wigner energy forms for all even-even nuclei with ^4=80. The obvious 
discrepancy can be observed by using two forms for Wigner energy. The solid triangles and 
solid circles denote the results with form (1) and form (2), respectively. The value of the 
extracted symmetry-energy coefficients Osym is larger with form ( 2 ) than that with form ( 1 ), 
especially for the range of I close to zero, and the discrepancy decreases with increasing 
isospin asymmetry I. It is therefore necessary to determine the Wigner energy of nuclei for 
a better description of symmetry energy coefficient. 


III. THEORETICAL FRAMEWORK 


In semi-empirical mass formulas the Wigner energy usually is decomposed into two parts 
[33, 34] 


Ew{N, Z) = -W{A)\N -Z\- d{A)6N,znnp, 


( 8 ) 


where W{A) and d{A) are smooth functions of the nuclear mass number A. The hrst term 
on the right-hand side of Eq. ( 8 ) contributes to all A 7 ^ Z nuclei. The quantity equals 
1 for odd-odd nuclei and vanishes otherwise, and therefore the second term d{A) is nonzero 
only for N = Z odd-odd nuclei. The Wigner effect mainly stems from the hrst term in 
Eq. ( 8 ). By combining the hrst term in Eq. ( 8 ), the traditional symmetry energy term 
(A — ZY!A is replaced by T(T + x) term [35-38]. So the odd-odd symmetric nuclei are not 
considered in the following calculation. T = |T;j| = lAzZi jg isospin value of the nuclear 
ground state, and I = (A — Z)IA is the isospin asymmetry of a nucleus. Then one has the 
relation. 


1 ^ 

2 


(9) 


The symmetry energy term including the Wigner energy can be expressed as 


sym 


(A,T) = ^^T(T + x) = 


Aa. 


sym 


A 


T^ + 


4a, 


sym 


A 


Tx. 


( 10 ) 
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FIG. 2: (Color online) (a) Two forms Wigner energy as a function of isospin asymmetry /, and 
(b) the extracted symmetry-energy coefficients Osym by using two forms Wigner energy applied to 
all even-even nuclei with ^=80. The solid squares denote the result of excluded Wigner energy. 


Inserting Eq. (9) into Eq.(lO), we can obtain the symmetry energy per particle expression 
as the function of mass number A and isospin asymmetry I, 


e:yra(A, I) = 


-sym 


(A I) 


A 


_ * TZ _ TZ 

®sym'‘ CLsyraJ^ T 


2nsym^ I 


( 11 ) 


where e*y^(A, I) = esym{A, I) -b and = asym(l + jjp)- Osym is the symmetry energy 
coefficient be expressed as a function of mass number A and isospin asymmetry I. 2asymX 
denotes the Wigner energy coefficient, the value of x is not well determined from nuclear 
masses, x = 1 is associated with neutron-proton exchange interactions in SU(2) symmetry, 
while X = 4 corresponds to the full supermultiplet symmetry SU(4)[39]. The further discus¬ 
sion on the Wigner energy can be found in Ref. [40-43]. Here x as a parameter is introduced. 
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named the Wigner energy parameter. The x value has crucial effect on the symmetry en¬ 
ergy coefficient, since the symmetry energy is the summation of the traditional symmetry 
energy and the Wigner energy. The different x value denotes the different Wigner energy. 
Inserting Eq. (11) into Eq. (3) and replacing esym(^,/), the nuclear energy part 

per particle Eq.(3) becomes 


en{A,I) — eQ{A) + e* {A, I) 


2x 


— eo(^) + asym(^5-f)(l + ’ 


( 12 ) 


Inserting Eq.(7) into Eq. (12), and take the difference of e„(y4, 1) between two isobaric nuclei 
with same odd-even parity. Eq. (4) becomes 


Ae:« = en{A,I)-e4A,I,) 


= aSUP - in + aWjP - I}) + 


-syinv 


A 


(|/|-|/d) + 


2 C!-sym 3 ^ 


A 


.(|/|3-|/,|3). (13) 


where i=l, 2, 3, ..., n, aiyL = J — aggA The dependence of reference nuclei (A, Ji), 
{A,l 2 ), ... , and {A, In) can be canceled through the summation, and the average value 


Ae*yj^ of the difference of symmetry energy can be expressed as 


Ae 


sym 


“ r(^®sym + + ••• + 


n 

= e, 


''sym 
n 


n “ 


i=l 


-| /£. -j /6 

i=l i=l 

^20^ - - V |J,|) + f^!p^(|J|3 _ _ ^ |J^|3)^ 

fl < ^ 4 r) < ^ 


A 


i=l 


n 


(14) 


2=1 


when neglecting the microscopic shell corrections of nuclei, the result of Eq. (14) Ae|yj^ = 
e„(y4, /) — ^ X]r=i ^n{A, li) is obtained by the measured binding energy per nucleon of each 
series isobaric nuclei compiled in AME2012. By using the expression of the right-hand side 
in Eq. (14) and fitting Ae*yj^ from more than 2200 measured nuclear binding energies, we 
obtain the optimal values J = 30.25 ±0.10 MeV, Oss = 56.18 ± 1.25 MeV, aiyL = 8.33 ± 1.21 
MeV and x = 2.38 ± 0.12 with an rms deviation of 23.42 keV. 












FIG. 3: (Color online)Symmetry energy coefficients of nuclei as a function of (a) nuclear mass num¬ 
ber A and (b) of isospin asymmetry I. The solid squares and open circles denote the experimental 


data and the fitting results by Eq. (14) with the optimum parameters values J = 30.25 

MeV, Oss = 56.18 MeV, agym = 8.33 MeV and x = 2.38. 


■( 1 ) 


IV. RESULTS AND DISCUSSIONS 


In Fig. 3 (a), we show the extracted symmetry energy coefficients of nuclei as a function of 
nuclear mass number. The solid squares denote the extracted symmetry energy coefficients 
from the measured nuclear masses by using in Eq. (13). The open circles denote 

the htting results by Eq. (14) with the optimum parameters values. One can see that the 
experimental value of obtained in our approach by Eq. (13) shows some oscillations 

and fluctuations, which is probably caused by the shell effects and other nuclear structure 
effects. In Fig. 3 (b), we show the same data as in Fig. 3 (a), but as a function of isospin 
asymmetry I. Form Fig .1 and Fig. 3 (b), we can End that the extracted symmetry energy 
coefficients depend on the corresponding isospin asymmetry of nuclei, which decreases with 
increasing isospin asymmetry I for the same mass number A, the largest values located 
in the range of nearly symmetric nuclei. However the Wigner energy parameter x value 
influences every parameters in Eq. (14). Fig. 4 shows the coefficients J, Oss, Osym (in MeV) 
and a deviation (in keV) as a function of Wigner energy parameter x. From Fig. 4 we can 
see that the coefficients J (solid squares), Ogg (solid circles) and aiym (solid triangles) increase 
hrstly then decrease with increasing x values in the range from 0 to 12. The rms deviation 
(J (down triangles) decreases firstly then increases with increasing x values. The minimum 
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FIG. 4: (Color online) The volume symmetry coefficient J, surface symmetry coefficient Ogs, the 
coefficient Ogym of term (in MeV) and a deviation (in keV) as a function of Wigner energy 
parameter x. 

value of cr = 23.42 keV is corresponding to the set optimal parameters values. One may thus 
expect the coefficient x to he somewhere between 1 and 4. The volume symmetry coefficient 
J ~ 30 MeV is insensitive to the value x in the range 1 < x < 4. The surface symmetry 
coefficient Oss is sensitive to the value x in the range 1 < x < 4, whose value changes from 
38.72 MeV to 65.85 MeV. The coefficient aiyh is more sensitive dependence of the value x in 
the range 1 < x < 4, from -6.98 MeV to 16.56 MeV. So we draw a conclusion from the hgure 
that aiyL is not well determined from nuclear masses since x is ill-determined. For example, 
we change x value somewhat from 1.5 to 1.6, the value of aiyL changes from negative to 
positive. So the sign (positive or negative) of UsyL is sensitively dependent on the value of 

X. 

The contributions of symmetry energy and Wigner energy are also studied. As an exam¬ 
ple, the contribution of per term is shown in Fig. 5, where the asymmetric nucleus Ji = 0.07 
is selected as reference nucleus, since it is the minimum value of known nuclei in mass table 
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FIG. 5: (Color online) (a) The values of Osym term (thin curve), Wigner term (thick curve) and 
flsym term (dashed curve) in Eq. (13) as a function of I, and (b) the contribution ratio of per term 
for A = 168. The parameters of J = 30.25 MeV, Oss = 56.18 MeV, oiyL = 8.33 MeV and x = 2.38 
are used. 

AME2012 for A = 168. From Fig. 5(a) one can see that the value of all three term increase 
with increasing isospin asymmetry I, when I < 0.39, the value of aiyL term aiyL(/^ — If) 
is less than that of Wigner term ^[aiyln(|/| — |/i|) + aiyln(|/|^ — |.^i|^)] in Fq. (13), when 
/ > 0.39 the value of OsyL term is larger than that of Wigner term. Fig. 5(b) shows the 
contribution ratio of per term, the ratio is calculated by the ratio of per term value to the 
summation of three term value. From Fig. 5(b) we can see the changing details of per 
term with increasing isospin asymmetry I. The average contribution ratio of four term are 
87.92%, 8.27% and 3.81% for aiyL term aiyL(/^ — If), Wigner term and aiyL term in the 
range of / = 0.07 — 0.5, respectively. With the increasing of isospin asymmetry /, the asym 
term is the major contributor, which increases hrstly and reaches a maximum at / = 0.27, 
and then decreases with increasing isospin asymmetry I. The Wigner term decreases and 
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the aiyin term increases with increasing isospin asymmetry I. The contribntion ratio of asyh 
term is less than that of Wigner term in the range of / = 0.07 — 0.39 and larger than that 
when / > 0.39. 

V. SUMMARY 

In snmmary, we have proposed a method to extract the symmetry energy coefficient (in¬ 
cluding the coefficient aiyL of term) from the differences of available experimental binding 
energies of isobaric nuclei. The advantage of this approach is that one can efficiently remove 
the volume, surface and pairing energies in the process. It is found that the extracting exper¬ 
imental symmetry energy e*y^(A, I) should be the summation of the traditional symmetry 
energy esym(^, I) and the Wigner energy ew(A, I). And a*yjjj(A, I) decreases with increasing 
of isospin asymmetry I, which is mainly caused by the Wigner energy effect. Through the 
polynomial ht to the result of Ae*yj^ by the right-hand side expression of Eq. (14), we have 
obtained the optimum parameters values J = 30.25 ± 0.10 MeV, Oss = 56.18 ± 1.25 MeV, 
dsym = 8.33 ± 1.21 MeV and the Wigner parameter x = 2.38 ± 0.12. We also hnd that the 
volume symmetry coefficient J ~ 30 MeV is insensitive to the value x, while the surface 
symmetry coefficient Ogs and the coefficient oiyL are very sensitive dependence of the value 
X in the range 1 < a; < 4, especially for aiyL, whose value maybe change from negative to 
positive since the change x value somewhat in the range 1 to 4. The contribution of the 
wigner energy term decreases and the contribution of asyL term increases with increasing of 
isospin asymmetry /. For very neutron-rich nuclei, aiyL term will play an important role 
since its contribution is larger than that of Wigner energy term. 
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